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Abstract 

For type IIB supergravity with a running axio-dilaton, we construct bulk so- 
lutions which admit a cosmological background metric of Friedmann-Robertson- 
Walker type. These solutions include both a dark radiation term in the bulk as 
well as a four-dimensional (boundary) cosmological constant, while gravity at the 
boundary remains non-dynamical. We holographically calculate the stress-energy 
tensor, showing that it consists of two contributions: The first one, generated by 
the dark radiation term, leads to the thermal fluid of M = 4 SYM theory, while 
the second, the conformal anomaly, originates from the boundary cosmological 
constant. Conservation of the boundary stress tensor implies that the bound- 
ary cosmological constant is time-independent, such that there is no exchange 
between the two stress-tensor contributions. We then study (de) confinement 
by evaluating the Wilson loop in these backgrounds. While the dark radiation 
term favours deconfinement, a negative cosmological constant drives the system 
into a confined phase. When both contributions are present, we find an os- 
cillating universe with negative cosmological constant which undergoes periodic 
(de) confinement transitions as the scale of three space expands and re-contracts. 
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1 Introduction 



Gauge/gravity duality pQ has proved to be extremely successful in describing strongly 
coupled systems. This applies in particular to confining theories which can be modelled 
for instance using non-trivial dilaton flows (2J EJ H]. In these models, the Wilson 
loop displays an area law. Moreover, gauge/gravity duality has also proved useful in 
describing deconfined finite temperature field theories which are naturally assumed to 
be dual to asymptotically AdS black holes. 

Both confinement and horizon formation also arise in quantum field theories on 
curved space backgrounds, in Anti-de Sitter and de Sitter geometries, respectively. In 
the gauge/gravity duality context, this has been investigated for instance in [5j[6j[71|8l 
[TUl [HI [12] by considering a boundary cosmological constant A in the four-dimensional 
boundary quantum field theory. Holographic studies of strongly coupled quantum field 
theories in curved backgrounds are however not only interesting in their own right (e.g. 
in order to verify the properties of particle production phenomena such as the Unruh 
effect at strong coupling), but also from the point of view of AdS/CFT dualities for 
time-dependent backgrounds. In particular, standard cosmological evolution in the 
presence of a cosmological constant can yield de Sitter or Anti-de Sitter geometries. 
Thus from studying gauge theories in these backgrounds we expect to learn about the 
properties of matter in the early universe (e.g. during inflation). 

With this situation in mind, in this paper we consider gravity duals of field theories 
on cosmological backgrounds where in the dual gravitational description a bulk radia- 
tion term is present in addition to a boundary cosmological constant. This term has 
first been considered in brane world models in [T3(, [14j HSl [16]. Due to its schematic 
form C/a 4 , with a being the scale factor, it corresponds to a relativistic radiation con- 
tribution to the energy density. We discuss the interplay between this radiation and 
the boundary cosmological constant in the boundary energy momentum tensor, as well 
as their effects on the temporal Wilson loop. We find that the combined effect of the 
dark radiation term and the boundary cosmological constant introduces an effective 
dynamics into the dual field theory, triggering (de) confinement transitions for the Wil- 
son loop. For vanishing boundary cosmological constant and flat horizon topology, 
the dark radiation term gives just the Stefan-Boltzmann contribution p ~ T 4 to the 
boundary energy density. In the other cases the relation to the temperature is more 
involved due to the time-dependence of the background geometry, as we discuss. 

As a further ingredient, we consider a running axio-dilaton similarly to the model 
of Liu and Tseytlin [I]. The axio-dilaton introduces a finite gluon condensate which 
on flat space leads to confinement. 

Our main result are explicit evaluations for Wilson loops in the field theories dual to 
the gravity solution with dark radiation term for the three cases of positive, vanishing 
and negative boundary cosmological constant. The general intuition arising from the 
static quark-antiquark potential is that the dark radiation term always drives the sys- 
tem into a deconfined phase (with the Wilson loop displaying a perimeter law), since 
it acts similarly to a temperature in flat space. We find an interesting pattern for the 
(de)confining behaviour of the Wilson loop, depending on the sign of the boundary 
cosmological constant: 
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1. For positive cosmological constant, the theory is always in a deconfined state, 
even for vanishing dark radiation term, which is in accordance with the general 
expectation that de Sitter-like expansion tends to destabilise bound states. 

2. For vanishing cosmological constant our bulk metrics are diffeomorphism equiva- 
lent to topological AdS black holes [I7J[T8]. In this case, the running axio-dilaton 
as in [4] is crucial for determining the confinement properties. The Wilson loop 
shows confinement if the dark radiation constant is vanishing and deconfinement 
otherwise, both for flat and hyperbolic horizons. As discussed further below and 
in sec. the non-trivial dilaton flow is essential for the Wilson loop confinement 
in this case, as in the absence of the running dilaton the quark- ant iquark poten- 
tial would be screened for a gravity dual with hyperbolic topological black hole 
[T5] at all temperatures. 

3. For negative cosmological constant, we find an interesting (de) confinement transi- 
tion which occurs due to the competition between the deconfining dark radiation 
term and the confining nature of Anti-de Sitter-like contraction: For small scale 
factors the Wilson loop is deconfined, while for large scale factors it is confined. 
Intuitively this can be thought of as the Wilson loop probing the (holographically 
defined) field theory vacuum whose energy-momentum VEV now has two com- 
ponents - the conformal anomaly component due to the boundary cosmological 
constant, and the dark radiation component which in essence behaves like ther- 
mal relativistic radiation, getting diluted by the scale factor as a~ 4 (t)|]] Hence 
close to the singularity the dark radiation component dominates, driving the sys- 
tem to deconfinement, while away from it the confining nature of the negative 
cosmological constant dominates. 

In this work we mainly focus, for simplicity, on the Wilson loop as a measure of quark- 
antiquark (de) confinement. On curved space-times, other measures of confinement 
such as the density of states or the mass gap criterion do not necessarily coincide with 
the Wilson loop criterion: For example, in j22j [23] it was argued that the Wilson loop 
confines on AdS spaces at any temperature, while [12] showed that Neumann boundary 
conditions on the boundary of AdS space allow for a large N deconfinement transition at 
finite temperature, with the deconfined phase being characterised by a 0(N 2 ) density 
of states at low energies. We hence leave a thorough investigation of the subtleties 
involved in relating these different criteria for confinement for future work, and rather 
focus on the Wilson loop as a criterion to characterize our holographic backgrounds. 
The reader should also consult sec. [5] for a more in-depth discussion of the case of 
vanishing cosmological constant and hyperbolic horizon: In this case the Wilson loop 
and the density of states measure indeed do not agree, due to the presence of the gluon 

^his interpretation of the dark radiation constant has also been given in the gauge/gravity context 
in [2Qj [TH Q2] ■ Furthermore, the works [HI [18] write the AdS-Schwarzschild black hole (without a 
running dilaton) in a cosmological foliation, which is possible in our construction as well. However, 
our setup is more involved due to the presence of the boundary cosmological constant and the running 
dilaton (see sec. [5]). Finally, a related study [21] for k = and without dark radiation found an interplay 
between boundary cosmological constant and the tension of an IR brane sourcing the geometry. 
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condensate (i.e. the running dilaton) and since the thermodynamical contributions of 
the axio-dilaton cancel each other in the Liu-Tseytlin ansatz. 

On the technical side, we decouple the axio-dilaton dynamics from the five-dimensio- 
nal metric in the same way as in [4]. We then solve the Einstein equations of five- 
dimensional Einstein-Hilbert gravity with a cosmological Ansatz for the metric already 
used in [TJ] in the context of brane-world cosmology. The Friedmann equation arises 
from the constraint equation of the bulk Einstein equations, and we include the dark 
radiation term into our analysis. By imposing the usual Dirichlet boundary condi- 
tions of holography, the four- dimensional boundary gravity remains non-dynamical, 
as there is no four- dimensional Einstein-Hilbert action on the boundary. Besides the 
dark radiation term we also allow for a boundary cosmological constant in our Fried- 
mann equation. We find that requiring boundary diffeomorphism invariance leads to 
a time-independent boundary cosmological constant. 

The remainder of this paper is structured as follows: In the next section, the holo- 
graphic background is given and we discuss how the boundary cosmological constant 
arises. In section 3, the holographic interpretation of the dark radiation term in our 
approach is illuminated using the specific example of vanishing boundary cosmolog- 
ical constant. Section 4 discusses the solution for finite boundary cosmological con- 
stant, and how boundary diffeomorphism invariance (i.e. conservation of the boundary 
stress-energy tensor) forces the cosmological constant to be actually time- independent. 
Section 5 then derives the main result of this paper: The Wilson loop expectation val- 
ues are calculated and their (de) confinement properties are classified. Summary and 
discussions are given in the final section 6. 

2 The Background Geometry 

In this section we first review the reduction of ten-dimensional type IIB supergravity 
to a five-dimensional dilaton gravity by a Freund-Rubin Ansatz which allows for a 
nontrivial axio-dilaton. This ansatz, first employed in [3j HJ 121], links the axion with 
the dilaton in a way which allows to describe 1/4 supersymmetric D3-D(-l) solutions. 
Supersymmetry can then be broken by introducing finite temperature. In 12.21 we then 
solve the five- dimensional Einstein equations with a time-dependent ansatz for the 
metric along the lines of [15] and find holographic backgrounds describing a cosmolog- 
ical evolution at the boundary. In this course we identify the boundary cosmological 
constant, driving the cosmological evolution of the boundary metric, when solving the 
constraint equations in the bulk. 

2.1 Five-dimensional Dilaton Gravity from IIB Supergravity 

We start from the ten-dimensional type IIB supergravity retaining the dilaton $, axion 
X and selfdual five form field strength F( 5 ), 

S = ±J (r - i(9*) 2 + ie 2 *(9x) 2 - j^ffo) , (1) 
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where other fields are consistently set to zero, and x is Wick rotated [23]. Under the 
Freund- Rubin ansatz for F( 5 ), F Ml ... M6 = — vA/2 e^...^ [3j S], and for the lOd metric 
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taken as M 5 x S 5 



^ s io = 9mn(Ix m dx N + g a bdx a dx b , M, JV = 0, . . . , 4 , a, b = 5, . . . , 9 , 
the equations of motion of the non-compact five- dimensional part M5 become □ 

Rmn = ]z (d M <$>d N $ - e 2 *d M xd N x) - ^9mn (3) 
1 --9m (V^9 MN d N ^) = -e™g MN d MX d N x , (4) 



d M (V=ge 2 *g MN d NX ) = (5) 

These equations have a supersymmetric solution when the following ansatz is imposed 
for the axion x [231 E], 

X = -e - * + Xo ■ (6) 
In this case, using the ansatz (jBJ) in (J3j) (J5j) gives rise to the two equations 

Rmn = —AgMN (7) 

and 

9m (V=99 MN d N e*) = , (8) 

where @ and (jSJ) now may be shown to coincide using (JSJ). The latter set of equations 
is also useful for finding finite temperature solutions in which supersymmetry is broken. 

2.2 Solution with Dark Radiation 

We examine here time-dependent solutions which include a "dark radiation" term 
[T5l IT6] (also known as 'mirage energy density' [20J). To find this term, we change the 
radial coordinate r to y, where r/R— fxr — e m and fi — l/R — vA/2, and we consider 
the following Einstein frame metric, 

ds| = —n 2 (t,y)dt 2 + a{t,y) 2 ^ i jdx l dx i + dy 2 , — . . . , 3. (9) 

In this metric, we obtain from the Einstein equation for the tt and yy components [15] 

a V k A fa'V C , . 

— +- = -T+ - 10 

na / a" 1 4 V a / a 4 



2 The five-dimensional part M5 of the solution is obtained by solving the following reduced Einstein 
frame 5d action, 

S=^Jd 5 x^—g(R + 2,K-\{d<i>f+ 1 -e^{dx) 2 ^ ■ (2) 

The opposite sign of the kinetic term of \ is due to the fact that the Euclidean version is considered 
here [Ml. 
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where a = da/dt and a' = da/dy. Note that this is a first order equation, integrated 
from the second order Einstein equations (see [15] for their explicit form). It then 
turns out that without any additional matter in the bulk, the integration constant C 
must be a constant with respect to both y and t in order to satisfy both the tt and 
yy components of Einsteins equations. This constant C appears in the equation (TTUT) 
in the form ^, which is usually referred to as "dark radiation" term, since it behaves 
exactly as a component of relativistic radiation which, in the context of braneworld 
models, leaks from the bulk into the UV brane |16j . 

It also needs to be checked whether the Bianchi identities and the ij and ty compo- 
nents of Einsteins equations are satisfied with the above Ansatz. As shown in [15] , the 
first two are satisfied upon the use of eq. ( ITOi) , while the latter relates the free function 
n(t,y) to a(t,y) up to a time-dependent integration constant, 



n' a a' 



o = ----. (11) 

n a a 

This last equation is solved by setting the following ansatz [151 116] - 

n(t, y) = , a = a {t)A{t, y) . (12) 

ao[t) 

Then the equation for A(t, y) is obtained from (ITU]) as 

"oV k A 2 2 . C 



-i A+(A)+ ^> (13) 

where A' = dA/dy. Looking at eq. (fT3|) . we recognize its left-hand side as part of 
the Friedmann equation from standard cosmology. More precisely, it is the part of 
Friedmann's equation without the cosmological constant term. In particular, since 
the left-hand side of eq. ( I13p is a function of time only, the right-hand side of (TT3"j) 
must also be only a function of time, i.e. independent of the radial coordinate y. 
The right-hand side of ffT3]) thus effectively acts as a time-dependent vacuum energy 
"source term" for the cosmological evolution at the boundary, described by the left- 
hand side. Introducing a time-dependent boundary cosmological "constant" A(t), we 
can thus separate eq. (fT3"j) into its left and right hand sides, yielding two independent 
equations. This procedure is similar to separation of variables when solving differential 
equations: For general time-dependent X(t) the above replacement of eq. ( JT3"j) by the 
two equations ( JT4")) and (|T5|) does not affect the solution space, since every solution of 
(|T4|) and ( IT5|) will be a solution of ( !T3|) . and vice versa. 

Doing so, the left hand side of ( !T3|) then becomes the four-dimensional Friedmann 
equation with a four- dimensional boundary cosmological term A(t), 

^y+4=A(t), (14) 



where k = ±1 or 0. From standard cosmology, (1141) only yields universes with spherical 
(k = +1) topology for A > 0, while for k = the allowed choices are sgnA = 0,+l, 



6 



and for negative spatial curvature k = —1 even a spatially homogeneous and isotropic 
universe of constant negative curvature is allowed, i.e. sgnA = — 1,0, +1 are possible 
choices. 

For any A(t), A(t,y) can then be solved for by the following first-order differential 
equation in the variable y, 

A(() = _^ + (^ + JL, (15) 

using the solution a (t) of (|T4|) . In the above treatment of eq. (TL3I) we introduced an 
a priori time- dependent function X(t). In an evolving universe, a time-dependent cos- 
mological constant however would have to be sourced by additional energy-momentum 
sources at the boundary or in the bulk, which generate the relevant piece in the energy- 
momentum tensor that ensures energy-momentum conservation. Since in the holo- 
graphic context with standard Dirichlet boundary condition, gravity at the boundary 
is not dynamical (i.e. the background metric for the dual field theory is a fixed back- 
ground field), no boundary matter source can influence the boundary metric. The 
holographic energy-momentum tensor itself thus has to be conserved. We will cal- 
culate the holographic energy-momentum tensor in sec. HI but quote the result here 
already and argue that stress-energy conservation forces the cosmological constant to 
actually be time- independent. 

The general solution of eq. ( TT5|) . which will be analysed in more detail in sec. HI is 



^(U-^V + lg^l • <i<0 




Using standard holographic techniques (for details see sec. I4.4.2p . the vacuum expec- 
tation value of the boundary stress-energy tensor is found to be of perfect fluid form, 

4i? 3 

(T^u) = diag(-p,p,p,p), a = ^ , (17) 

16nG y N ' 

( C \(tf\ ( C 3A(t) s 

P = 3a 77^71777 + ~^7T ' P = a 



4R 2 a*(t) 16 ) ' * \4R 2 a*(t) 16 

If we now impose the holographic stress-energy tensor to be conserved, V M (T M j,) = 0, 
we actually require a continuity equation for pressure and energy density, 

= p + 3H(p + p) . (19) 

In order to satisfy this continuity equation, the boundary cosmological constant A(t) 
then has to be a constant, 

A(t) = 0. (20) 

Physically speaking, requiring the stress-energy tensor to be conserved amounts to re- 
quiring the holographically defined generating functional to be invariant under bound- 
ary diffeomorphisms. The conservation equation V ^(T^) = then is the one-point 
function diffeomorphism Ward-Takahashi identity. In holographic renormalisation this 
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is a natural outcome since both the regularised on-shell action as well as the coun- 
terterms are constructed in a manifestly boundary diffeomorphism invariant way. A 
possible non-conservation of the bulk part of (T^) can then only be cancelled by addi- 
tional boundary terms which change the chosen boundary conditions from Dirichlet to 
Neumann or mixed ones [25] , hence inducing additional dynamical degrees of freedom 
into the boundary theory. We thus conclude that to constitute a physically meaning- 
ful holographic background, the boundary cosmological constant must be an actual 
constant in time. One should note that restricting the solution space of eqs. f fT3|) or 
(jl4l) . (115p does not influence the argument given above concerning the equivalence of 
the (restricted) solution space of both sets of equations. 

We will show in the following (in particular in section |5]) that the solution A(t,y) 
of this equation encodes important dynamical properties of the gauge theory in a 
Friedmann- Robert son- Walker universe. As will be shown in section [5j the behaviour 
of Wilson loop expectation values as calculated from a minimal string world sheet, and 
hence the (de) confinement properties of the vacuum show an interesting competition 
between the dark radiation constant C > 0, and the boundary cosmological constant 
A. 

Finally, we would like to comment on the physical situation concerning the dark 
radiation term and the boundary cosmological constant in brane world models [26| [27] , 
which is slightly different. In these models, due to the fact that the UV brane sits 
at a finite cutoff and due to the chosen boundary conditions, an energy exchange 
between the bulk and the brane is possible. In particular, the value of A is always 
tuned by the bulk metric and the five- dimensional cosmological constant A, and the 
dark radiation term is considered to be an energy flux between brane and bulk. 
In the holographic setup we consider here, due to the standard Dirichlet boundary 
conditions chosen in holography, no bulk-boundary energy exchange is possible, and 
the boundary cosmological constant can be freely tuned. The holographic setup is thus 
less rigid compared to the brane-world models. 

3 Holographic Interpretation of Dark Radiation 

Above we have shown how to obtain a consistent background solution describing a 
boundary metric undergoing cosmological evolution under the influence of both a 
boundary cosmological constant and a dark radiation term in the bulk. Here we con- 
centrate on the holographic interpretation of the dark radiation term for the simplest 
case, i.e. the case of vanishing boundary cosmological constant. We find that the dark 
radiation term introduces a temperature for the boundary = 4 field theory. 

3.1 Solution for Vanishing Boundary Cosmological Constant 

For the case of vanishing boundary cosmological constant, a solution of fl!4p is given 
by A = 0, k = and ao(t) = 1, and A(t,y) = A{y) is obtained by solving (|T5|) with 
A = 0. This gives 

A = e m (1 + c e-^ y ) l/2 , (21) 
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where Co = C/(4/i 2 aQ) = C/(4 / u 2 ) since ao = 1. From Eq. (fT2|) . we obtain 

1 2C 



A X 



Then, using r/R = e^ y the full Einstein metric is given by 



(22) 



ds 2 = i-n 2 dt 2 + A^dx') 2 ) + %dr 2 + R 2 d£ll , (23) 

For the dilaton (jSJ) we find, using the above metric, 

c * I, g log l+Mg/!j! f2 5) 
6 " i+ 25 0J R 4 g l-5o(W 1 j 

Here the integral constant g corresponds to the gauge condensate (trF 2 ), and the 
boundary condition e* — > 1 for r — > oo is imposed. For c = 0, this solution reduces 
to the supersymmetric one used in |28j, and e* diverges at r = if c > 0. 



3.2 Holographic Interpretation of Dark Radiation 

In order to give an interpretation to dark radiation constant C, we rewrite the solution 
eq. ( 1231) as a planar AdS-Schwarzschild black hole. The five- dimensional part of the 
metric in the Einstein frame can be brought into that form, 

d*h = ^(-/(^ 2 + (^) 2 ) + | M , m = l-£, (26) 



by the coordinate redefinition de 



1 + ^r = (CR 6 )^. (27) 

Thus the dark radiation constant C has to be positive, and sets the horizon radius of the 
AdS-Schwarzschild black hole. The dark radiation constant is nothing but the mass of 
the AdS-Schwarzschild black hole and, applying standard holographic renormalisation 
[29| |30| 131] we find the dual stress-energy tensor to be of perfect fluid form 

<T$> = -^diag(3, 1,1,1). (28) 
16nG N 

We thus conclude that in holography the dark radiation constant defines a temperature 
for the fields of the dual field theory. For a non-expanding cosmology (k = A = as 
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in this case) this directly leads to a field theory (in this case M = 4 with gluon and 
instanton condensate) at finite (Hawking) temperature 



-Ho 



(45o) 



IT R 



(29) 



Using this temperature and the energy momentum tensor ( 13.21) . we can in particular 
confirm the Stefan-Boltzmann law for the energy densitjip = -<T° ), 

4i? 3 ( ~cp \ SN 2 2 4 

l 3 ^h^ T - (30) 

where we used = 8rc 3 a' 4 g' 2 / R 5 and i? 4 = ArcNa' 2 g s . This expression reproduces 
the known results of (221 E3]- In brane- world models, the dark radiation term has 
been interpreted as the radiation of the bulk gravitons which transfer the energy of the 
fields in the brane to the bulk. In [3H [16] it was noted that the dark radiation constant 
corresponds to the mass of the bulk AdS-Schwarzschild black hole. In their contexts, 
gravity is dynamical on the UV brane, and the dark radiation term appears in the 
Friedmann equation on the brane. As noted already in section [21 the holographic setup 
considered in this work is different due to the Dirichlet boundary conditions imposed. 
Here, gravity is not dynamical at the boundary of space-time. Instead, in our case the 
dark radiation term is dual to the energy density of the M = 4 U(N) SYM fields in 
a thermal state, as evident from the Stefan-Boltzmann law fj30|) . The dark radiation 
constant C, which appears as an integral constant when solving Einsteins equations 
[To] , sets the temperature of the dual field theory. To the best of our knowledge 
such a holographic interpretation of the bulk radiation term has not yet been given 
in the literature before. This interpretation will qualitatively also hold in the time- 
dependent cosmologies considered in sections H] and [5j We find that in all cases the 
dark radiation constant contributes in a thermal manner to the holographic stress- 
energy tensor of the system, with a time-dependent prefactor a (t)~ 4 associated with 
the dilution of relativistic radiation due to expansion or contraction of the (boundary) 
universe. On the other hand, the boundary cosmological constant yields a conformal 
anomaly contribution to the stress-energy tensor. We will see that both contributions 
can compete, giving rise to interesting dynamics. 



4 Holography for Boundary (A)dS4 Space-Times 

Above we saw that for vanishing boundary cosmological constant, the dark radiation 
constant corresponds to a temperature for the M = 4 fields. In this section we treat 
the case of finite boundary cosmological constant, and discuss in particular the bound- 
ary stress-energy tensor. We show that the dark radiation term induces a relativistic 
radiation contribution to the boundary stress-energy tensor, varying in time with the 
well-known a (t)~ A dependence during cosmological expansion. Furthermore, we find 
that stress-energy conservation in the boundary theory forces the boundary cosmolog- 
ical constant to be time- independent. 

3 The background metric in this special case is just the flat Minkowski metric ry Miy = diag(— 1, 1, 1, 1). 
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4.1 Solution for Finite Boundary Cosmological Constant 

A solution of (fT5|) for finite \(t) is 



A 



-j'tj 



4/i 2 



2 \ V2 



(31) 



where Co = C/(4/i 2 ao(t) 4 ). Here we have chosen asymptotic boundary conditions 



oo 



r/R, 



(32) 



where fi = 1/R, i.e. we require the asymptotic form of the metric to be AdSs. We find 
that A has time-dependence through ao(t) in cq and also A(£). This point is important 
to determine the structure of the metric below. 
From Eq. (fl2|) . we obtain 



3 2/iJ/ 



n 



A(t,y) 



1 - 



ACt) 



Co(*)e" 



(33) 



With r/R — e w , the full Einstein frame metric is then given by 

.2 ^2 



ds 2 



(-n 2 dt 2 + A^alitWixMdx') 2 ) + -^dr 2 + R 2 dnl , 
it 2 H 



(34) 



where 



.4 




n 



1 - 



l/i- 



(f) 2 ) -5«(f) 4 



1 - 



4^ 



(?) + *>(?) 



(35) 

The above metric has no naked singularities for time-independent A, as we checked 
by calculating R, R^R^ and the Kretschmann scalar R lu , p<T R. llv,xr jQ We will see in 
sec. 14.41 that A also needs to be time-independent in order to ensure boundary energy- 
momentum conservation. Note that we use a coordinate system in which the constant 
curvature three-space has the metric 



dQ 2 



dx 2 



fc£2/4) 2 ' 



(36) 



These are simply the standard spherical coordinates on the isotropic and homogenous 
three-space, with a conformal factor. 



4 Recently it was noted in a similar but not identical construction [3 5) that naked singularities 
might appear when deviating from pure dS expansion. Their singularity so far cannot be shielded by 
a horizon. In contrast, the Einstein frame curvature singularities in the backgrounds considered here 
are always behind the horizon g tt = 0, and coincide with the cosmological singularities ao(t) = 0. 
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4.2 Almost Constant Scale Factor and Adiabatic Expansion 

Quantum fields in an expanding space usually are not in thermal equilibrium, not even 
locally, unless the expansion rate is slow compared to the equilibration time of the 
system. This should be the case for very small but nonzero boundary cosmological 
constant A, in which case the scale factor a Q (t) would still be changing with time, but 
with a very slow rate. In other words, the Hubble rate H = ao/ao is small. In this case 
we can still make statements about the "slowly varying" temperature of the system, 
corresponding to the adiabatic regime]^ We find the horizon as the zero of the gu 
metric coefficient in (Ell), which is at 




r H = Rj~ c ^ + — (37) 



for A > — (A/i 2 )c^/ 2 = A c . If the scale factor a is slowly changing, it is possible to 
approximately satisfy (|T4|) with a time-independent ao, by taking k = 1 for A > and 
k = —1 for A < in the Friedmann equation (jT4j) . 

ao^l/IAI 1 / 2 , >y(x)=(l + k^ . (38) 
In this case, <9 T r# ph 0, and from the near-horizon geometry 



^.ggj 2 ^^ R 2 de 2 + . . . _ (39) 

a (slowly varying) Hawking temperature can be found for A > — (4/x 2 )\/co, reading 

VX/W) + (45 )V 2 
Th = • (40) 

We thus find that negative (positive) A decreases (increases) the effective temperature 
for the dual field theory. Furthermore we observe that the regime A < — (4/i 2 )\/co is 
special: Formally, the Hawking temperature calculation does not apply to that case 
even if dark radiation is present, since g# has no real zero any more, i.e. there is no 
horizon. We will see in sec. |5]that in this regime the Wilson loop shows a confining area 
law behaviour. The situation is thus similar to the Sakai-Sugimoto model [36], where 
the gravity dual of the confined phase is a cigar-shaped geometry which smoothly caps 
off instead of admitting a black hole horizon. 

4.3 Dilaton Solution 

In addition to the metric considered above, the other important field in our system 
is the running dilaton, whose exponential is related to the gauge coupling in the dual 



°The system evolves adiabatically, starting from to, roughly for a time span <J \\\\t — to| < C 1 
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field theory. The solution to the dilaton equation of motion reads 



Q 



2c 1 + 



A 2 



16coai 4 



log 



1 + d (R/r) 4 + (\R/(4y 2 r)) 2 ((R/r) 2 - 8/i 2 /A) 
1 - ~c Q (R/ry + (Ai?/(4/i 2 r)) 2 (( J R/r) 2 - 8/z 2 /A) 



+ 



A 1 — i 

(tan -1 /3 + tanh" 1 — tt)^+7. 



where q and 7 are the integration constants and 

(r/R) 2 - A/(V) 



(41) 



xV2 



(42) 




(1.5 



1.0 



1.5 



2.5 



Fig. 1: Plots of vs r for (a) A = -1 - 4(i 2 c /2 , (b) A = -^ 2 c /2 and (c) A = 
1— 4/i 2 Cq / ' 2 . Cases (a) and (c) are taken as examples for A < — An 2 c\j 2 and A > — 4/i 2 Cg / ' 2 , 
respectively. Other parameters are set as = R — 1, q = 2 and c = 0.1. In the case 

of (c), diverges at the horizon r# 



^ /2 + A 



which is 0.5 in this case. Note 



-0 1 4/^ 

that in general Co is explicitly time-dependent, these curves represent snapshots of the 
dilaton solution at constant time. 



We notice the following points for the above solution ( l4Tj) : 

1. The above expression (141 p seems to be complex due to the factor in the 
second line of f )4ip . However, this is necessary to cancel the imaginary part of 
tanh -1 (3 which has a constant imaginary part in/2 for (3 > 1. The condition of 
(3 > 1 is realized for r > r# (all r > 0) in the case of A > -A^cJ 2 (A < -4/i 2 Co /2 ), 
and hence we only give the solution of e* in this regime, and explicitly display 
the factor ^vr. As a result, the above expression (|4"T|) is real. 

2. The factor arctanh(/3) diverges for — > 1, which is realized for r — >■ r# = 

-^l/ <</ 2 + 4^2"- The same logarithmic divergence comes from the first logarithmic 
term in the equation (l4"Ij) . This divergence can be seen in case (c) of Fig. [TJ 
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3. For the case of A < — 4^i 2 Cq //2 , the solutions extend to r = and there is no 
divergence at any point in radial direction. The value at r = is given by 

$(o) . 8g// A „ ^ + 16/i 4 c 



A 2 + 1Q/j,% V. A 2 - 16/i 4 c 
A 



2^V 



(tan -1 Pq + tanh -1 /3 — vr) , (43) 



where /?o = A/(4cg // 2 ). The important point is that e*^ ^ is finite. We plot its 
numerical value in Fig. [2] for an appropriate parameter set as a function of |A|. 

It is interesting to note that the asymptotic value of e*^ at A = ±oo is given 
by 7. Thus, in the limit of asymptotically large positive or negative cosmological 
constants, there is no running of the coupling due to the gluon condensate, but 
only due to the conformal anomaly induced by the background. It would be 
interesting to further investigate this fact from a field-theoretic point of view. 

4. Further, if we set 7 = 1, we obtain the following asymptotic form 

~ 1 + q/r 4 + ■■■ (44) 

as r — > 00. This is the standard AdS/CFT expansion for a scalar dual to a 
A = 4 operator, which in this case is the gluon condensate TrF 2 . The integration 
constant 7 corresponds to the nonnormalisable mode, while q encodes the vacuum 
expectation value (TrF 2 ). 

Thus, while the ultraviolet behaviour of e* does not depend on A, the behaviour 
near the infrared region is very sensitive to the boundary cosmological constant. For 
A > — (4/! 2 )a/co, e* diverges at the horizon of the black hole configuration. On the 
other hand, for A < — (4/i 2 )y / co, e* approaches a constant at r = 0, and d r e \ r= o = 0. 
In this case, then, the Yang-Mills coupling constant reaches at an IR fixed point for A < 
— (4fi 2 )y^CQ. We should however note that conformal invariance of the boundary theory 
is still broken due to the gravity contribution to the conformal anomaly. However, this 
will not affect on the renormalization group equation for the Yang- Mills partjj We 
thus would naively expect quark confinement for A > — (4/i 2 )-y/co due to the strong 
infrared coupling. However, in this case the Wilson loop calculation of section [5] shows 
that the quarks are not confined: The Wilson loop deconfines due to the presence of 



3 More exactly, the RG equation for the effective action T reads 



The right hand side vanishes since the Eulcr density e fJ,,yal3 e f " 7lS R^po-Rap^s is topological and vanishes 
when integrated over space-time, and the Weyl tensor Cp, vpa — since the cosmological backgrounds 
are conformally flat. Hence, if no operators except the gluon condensate are present and Pym in 
the IR, the theory approaches an IR hxpoint when the dilaton approaches a constant in the infrared. 
This is the case in the Liu-Tseytlin like backgrounds considered here, since [37] there the beta function 
vanishes in spite of the presence of a gluon condensate. 
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the horizon at gu = 0. We find confinement for A < — (4/i 2 )^/co instead, where the 
coupling constant is finite and not so large, but where the horizon is absent. We thus 
conclude that similarly to the situation in the Sakai-Sugimoto model [36], the main 
factor controlling the confinement dynamics in this setup is not the coupling constant 
(i.e. the running dilaton) but the presence of a horizon in the bulk. The dilaton 
running is, however, important at zero temperature (C = 0) and leads to Wilson loop 
confinement e.g. in the hyperbolic case k = —1, as will be discussed in sec. [5j 



—* 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ' I A| 

2 4 6 8 10 

Fig. 2: Plot of e^ ** vs |A|, where the parameters are set as in Fig. [H yielding AjjL 2 c l J 2 = 
1.265. We notice that is real for |A| > A/j, 2 c^ 2 = 1.265, which corresponds to case 
(a) in Fig. [H with no horizon present. For |A| < Ap 2 ^ 2 = 1.265, case (c) of Fig. [TJ the 
dilaton becomes complex in the region hidden behind the horizon. Moreover, — > 1 
for |A| — > oo. Note that since in general en is explicitly time-dependent, this curve 
should be understood at a given instance in time. 



4.4 Boundary Energy-Momentum Tensor and Boundary Dif- 
feomorphism Invariance 

4.4.1 The VEV of the Boundary Energy Momentum Tensor 

Next, we calculate the four- dimensional stress tensor from holography. The Fefferman- 
Graham expansion of the metric ( 1341) reads 



ds = — + , [i,u = Q,...,3, (45) 



dp 2 g^u^PjX^ 1 ) dx^ dx 
4p 2 p 

9?u(p, a^) = 9{o)nv + 9(2)^9 + P 2 (9(4)^ + h miiu log p + h 2{4)iJLV (\og p) 2 ) H (46) 

g(o)fiv = (0(0)00, 9(o)ij) = (-1, ao(£) 2 7ij) , 9(2)^ = ~ 2R 2 u? 9(f3)>J ' v ' ^ 



and 



48c -A 2 //i 2 16c + A 2 //i 2 
5(4)00 ~ 16^ ' 9{4)ij ~ 16^ 9{0)ij ' ^ 
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Then by using the general formula [29] 



( T M"> = lfaG N ~ \ 9 ^ V (( T ^( 2 )) 2 " Tr #(2)) " \ (9(2))^ + ^0(2)^1^/(2)) 



we find the holographic stress energy tensor 



(49) 



4/? 3 /— 3A 2 \ 

<T "" > = <f - } + ^ (V**-) • (50) 

4 f? 3 r 

« } > = TT^m j?&9m)> ( 51 ) 

where (T$/) is the "thermal" stress tensor contribution, i.e. the contribution which 
would be thermal for the M = 4 SYM fields if the universe was not expanding. The 
second term, which depends on A, comes from the loop corrections of the SYM fields in 
a curved space-time, and is the conformal anomaly contribution. The first term does 
not contribute to the conformal anomaly as in the usual finite temperature case. The 
conformal anomaly then reads 

q\2 

TO = -^ iV2 ' (52) 

where we used G$ = 87r 3 a /4 g 2 / 'R 5 and R 4 = AnNa' 2 g s . The conformal anomaly 
precisely matches the free field theory result which can be easily obtained using the 
general formulae given in [381 EH]. This is expected since the conformal anomaly is 
one- loop exact, and hence trivially interpolates from weak to strong coupling. 

4.4.2 Time Independence of the Boundary Cosmological Constant from 
Diffeomorphism Invariance 

The stress-energy tensor ( 15 Op has, in flat space (k = 0), the form of a perfect fluid 
stress-energy tensor. As such, it obeys the continuity equation 

p + 3H(p + p) =0, (53) 

where p, p and H represent the energy density, pressure and the Hubble constant 
H = d /a , and dot denotes the time derivative. From (J50]l we read off 

p = Sa \ — - H , p = a \ — - — S — , a = 7=r . 54 

A priori, from the way how we introduced the boundary cosmological constant X(t) 
in sec. 12. 2\ it has to be considered as being time-dependent. The continuity equation 
then requires 

A = , (55) 

where we used Co = C/(4p 2 ao 4 {t)). Continuity of energy density and pressure thus 
dictate the cosmological constant to be an actual constant in time. 
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The requirement that the boundary stress-energy tensor satisfies the continuity 
equation (|53|) is very natural from the point of view of diffeomorphism invariance of 
the boundary theory. Gauge/gravity duality provides in particular a way to holograph- 
ically calculate the generating functional of correlators in the dual field theory via the 
Gubser-Klebanov-Polyakov-Witten relation [JJ from the appropriately renormalised on- 
shell gravity action. Since the regularised on-shell action as well as the holographic 
counterterms are invariant under boundary diffeomorphisms by construction, 

the resulting generating functional respects this invariance in the absence of external 
sources. The boundary diffeomorphisms are the ones compatible with the Fefferman- 
Graham expansion of the metric, and hence can only depend on boundary coordinates 
x p . In particular, only the leading piece of the expansion, the boundary metric, trans- 
forms under boundary diffeomorphisms as gffl >->■ g$ + d^ u ), but all the subleading 

coefficients such as the stress-energy tensor VEV g$ will be invariant. The source- 
operator coupling J dPxg^^T^ then enforces the Ward identities 



upon transformation of the effective action by such a boundary diffeomorphism. This 
identity should hold for correlators involving the stress-energy tensor. In particular, 
the holographic stress-energy tensor should be conserved, 



For a perfect fluid T^ u = di&g(—p(t),p(t),p(t),p(t)), eq. (j5"7|) is equivalent to the con- 
tinuity equation fl53|) . The time-independence of A, eq. fl55|) . thus follows directly from 
the requirement of stress-energy conservation. 

In a more general setup one would however expect the boundary cosmological con- 
stant to change with time due to energy exchange between bulk and boundary. As 
noted before, this is not possible in the holographic setup, due to the imposed Dirich- 
let boundary conditions imposed at the boundary, and due to the fact that gravity 
decouples from the gauge theory in this the UV cutoff is taken to infinity [34]. 

This dictates that the energy-momentum contributions to the holographic stress-energy 
tensor coming from the bulk (the dark radiation part) as well as from the nontrivial 
boundary geometry (the boundary cosmological constant part) cannot mix with each 
other in a nontrivial (time-dependent) way. Hence, each of them is conserved by itself, 
yielding f )55|) . 

5 Dark Radiation, Boundary Cosmological Constant 
and Quark Confinement 

In this section we consider the combined effect of both the dark radiation term and 
the boundary cosmological constant on infinitely heavy quarks, i.e. test quarks, in 
the Super Yang-Mills theory, by holographically evaluating the static quark- ant iquark 
potential from a Wilson loop vacuum expectation value. One way to introduce (super- 
symmetric) quarks in the present context is through probe D7 branes [3D]- The test 



= V^(T^...) 



(56) 



= V M (T^}. 



(57) 
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quark-antiquark pair is then described by a string worldsheet ending on a prescribed 
space-time contour on the D7 brane. If the D7 brane corresponds to an infinitely mas- 
sive embedding, and barring special issues such as the presence of gauge field charge 
on the brane, the brane embedding then will coincide with the asymptotic boundary 
of our space-time, and we can consider a string worldsheet ending on a contour at this 
boundary. The static quark-antiquark potential is then calculated from the energy of 
the string, evaluated on a minimal surface with the boundary condition prescribed by 
the contour [41J. Usually, the string worldsheet then has two possible configurations: 

1. A pair of parallel strings, which stretch between the boundary and the horizon. 
This configuration describes a free quark-antiquark pair, and corresponds to a 
deconfined situation in which the Wilson loop shows a perimeter law. 

2. A U-shaped string whose two end-points are on the boundary, but which does 
not touch any black hole horizon or singularity in the bulk. If the g xx component 
of the string frame metric has a minimum, the string will be stuck there and the 
energy will depend linearly (for large separations) on the separation of the end- 
points, and show an area law for the Wilson loop. This configuration describes a 
confined quark-antiquark pair. 

These two types of configurations are seen to compete thermodynamically in the finite 
temperature gauge theory [11], as well as for the theory in dS4 [9], with the deconfined 
configuration being thermodynamically preferred in both cases. 



5.1 The Wilson Loop in Cosmological Evolution 

Following [H], we consider the Nambu-Goto string dynamics with the string world 
volume in (t, x) plane. The energy E of this state is then obtained as a function of the 
proper distance L between the quark and antiquark as follows Choosing a gauge 
X° = t = t and X 1 = x l = a for the world sheet coordinates (r, cr) of the Nambu-Goto 
action, the Nambu-Goto Lagrangian in the present background f[2"3"j) becomes 



J da e^ 2 n(r)\/r>i + (Ij* (A(r)a (t) 7 (x)) 2 , (5f 



1 f i <I>/2- / \ I ,o , ( r \ 



where only the radial coordinate r(x) is assumed to depend on x, and prime denotes 
the derivative with respect to x. The functions n, A are defined in f )34p and ( 13 5p . 

The energy of the string configuration, which is nothing but the static quark- 
antiquark potential, is obtained from fl58|) as 



1 /■._,,/. ( R 2 s 1 



where 



£ =- iNG = 2^ d °^t + \—A d * r] ' (59) 



a = a (t) I daj(a) = a (t) \ do - - ^ 2 , (60) 
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and 



,*/2 



\An\ 



(61) 



We should note that a measures the physical length (proper distance) in the present 
case, while a = x measures the distance in comoving coordinates, which do not change 
with expanding scale. We will consider the static quark- ant iquark potential as a func- 
tion of proper distance. 

The quark-antiquark potential fl59l) shows a scaling with distance if n s (r) has a 
finite minimum at some distance r = r* outside the horizon. This is the case for 
A < -VgJ /2 . In this case the dilaton in n s (r) varies very slowly and monotonically (see 
the discussion around figure [1]), and we can estimate the minimum of n s by neglecting 
the dilaton dependence, e*/ 2 ~ 1, in n s (r), i.e. taking 



r \ ■ 



r) M 



r y 
RJ 




c 



R 



(62) 



We then find the minimum of n s (r) at 



r* = R 



.V 



i/i 



- c 



(63) 



We see that the minimum is at a finite value of r* for A 2 > (4/i 2 ) 2 c > 0, since c"o > 
is necessary in order to have a positive temperature contribution of the Yang-Mills 
fields to the holographic energy-momentum tensor (I5TT) . There are thus two regimes, 
A > 4/i 2 \/co and A < —A^y/c^. In the former case, however, the dilaton e* diverges 
at a finite radius, and cannot be neglected any more in (16T|) . We are thus left with 
considering the case A < — 4/x 2 y 
is 



For A < — 4/x vcq, the value of n s at the minimum 



riAr 



2i? 




-£ >o. 



This is finite since we are considering the case of A < - 
since c"o > 0. Then the energy E is approximated as [9] 



-VcT. 



Note that n s (r* 



where 



2-rca' 



^O" max 

L = 2 I da 



(64) 



> 



(65) 



(66) 



is the proper distance between the string endpoints, and a m i n (cr max ) is the value at 
r min ( r max) of the string configuration [10J. The potential V — E thus grows linearly 



in proper distance as long as n s (r*) > 0, and the string tension is given by 



'qq 



n(r*) 
2-na' 



(67) 
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Fig. 3: Plots of n vs r for (a) A = —1 — 4 / u 2 c / , (b) A = — 4/i 2 c / and (c) A = 
1— 4/i 2 Cg //2 . Cases (a) and (c) are taken as examples for A < — 4fi 2 c^ 2 and A > — 4/i 2 c^ 2 , 
respectively. In agreement with ( 1631) . (a) and (b) show minima and hence confine. The 
parameters are taken to be = R = 1, k = —1, and cq = 0.2. For (c), there is a 
horizon at r = 0.5. Note that for the parameter values chosen in this plot, A < in 
all cases and hence k = — 1 is consistent when solving the Friedmann equation (fl4"|) . k 
only enters the dilaton and hence the Wilson loop via the value of Co- Furthermore, 
since in general c is explicitly time-dependent, these curves represent snapshots at 
constant times. 



On the other hand, for the case of A > — 4/i 2 c ( 1 / 2 we do not find such a finite 
minimum of n s (r). In this case n s (r^) = at the horizon r# defined by n(rn) = 
0. The exact behaviour of n s , including the dilaton dependence, is shown in the 
numerically obtained plot Fig. [3] for all cases of relevance. The numerical results 
support the approximations made above. In particular, there is no linear potential 
for A > — 4/i 2 CQ ,/2 . The quarks are deconfined in this case, which can qualitatively be 
understood as the effect of the cosmological constant not being sufficiently negative 
to overcome the thermal screening of the quark- ant iquark force. This is particularly 
interesting for negative A, in which case we find a possible deconfined phase even for 
AdS backgrounds if only the finite "temperature" screening of the quark- ant iquark 
potential, set in this case by the dark radiation constant C, is strong enough. This 
seems to be a novel phenomenon at strong coupling, since the Wilson loop in AdS 
spaces previously was expected to confine [221 123] due to the diverging gravitational 
potential in AdS space (for a more thorough discussion see section l5~2l ). 

We obtained the relation E(L), as shown in fig. HI and the tension T q g in the following 
way: Since the Lagrangian in f )58p does not explicitly depend on the coordinate a = x, 
we find the following quantity conserved under cr-shifts, 

e $/2 ; 1 (-YnA 2 (r)=H. (68) 

v/(r/i?) 4 A 2 (r) + (r') 2 \R' 

We can fix H at any point we like, so we fix it at r = r min . Then, choosing H = 
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e*/ 2 (£) 2 n{r)A{r)\ rmin , we get 



1 



/•rmax 

L = 2R 2 / cir — _ 

■Aw™ r 2 A(r)^/e*Wr 4 n(r) 2 i(r) 2 / {e^ r ™«hi in n(r min ) 2 A(r min y) - 1 

'"'"max ^iVV^O*)/ 2 

E = — I dr K ' (69) 

4 /i _ e*(^™)r^ n n(r miri )M(r mm ) 2 / (e $ Wr^(r)M(r)2) 



1 

7ra' 



Figure H] shows the exact dependence of the energy E on the distance L for the values 
q = (i.e. for constant dilaton) for A < — 4 / u 2 c ( 1 / 2 (curve A) and A > — 4/i 2 Cg' /2 (curve 
B). In the former case, we find the linear potential at large L as expected, and we 
find a typical screening behavior in the latter case, similar to the one seen in the finite 
temperature deconfinement phase. The qualitative behaviour of the Wilson loop for 
q 7^ are unchanged from the q = case. 



2.5 - 




0.5 l.O 1.5 2.0 

Fig. 4: Plots oiEvsL for (A) A = -l-4/i 2 gJ /2 and (B) A = +l-4/x 2 gJ /2 . Cases (A) and 
(B) are taken as examples for A < —A^, 2 c^ 2 and A > — 4/i 2 Cg / ' 2 , respectively. Case (A) 
shows deconfinement, while case (B) shows typical screening behaviour: In the latter 
case there exist two U-shaped configurations for small L, while for large separations 
the string breaks at the horizon and hence has vanishing energy (at 0(N 2 )). Other 
parameters are set as q = 0, // = 1, R = 1, k = —1, c = 1.0, r max = 3 and a' = 1 . 
Note that in general c is explicitly time-dependent, these curves represent snapshots 
of the dilaton solution at constant time. 
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5.2 Classification of (De)confining Behaviour in Cosmological 
Backgrounds 

In the above analysis we found confining behaviour for the Wilson loop if 



Thus for vanishing dark radiation constant C = (which would correspond to vanishing 
temperature in the static case) and negative boundary cosmological constant A < the 
system is always confined, in accordance with earlier results of J22J [12]. On the other 
hand, the Wilson loop cannot have an area law for 



This regime is deconfining. The right-hand side of these inequalities is generically 
time-dependent, and hence the inequalities will hold only in particular time intervals, 
in which the universe is larger /smaller than a critical value set by the boundary cos- 
mological constant. There are three cases: 

1. Positive A: In this case the system is in a deconfined phase. This could have 
been expected by the expanding nature of the de Sitter universe in this case: The 
expansion of space tends to destabilise bound states, leading to deconfinement 
even at zero dark radiation constant. 

2. Vanishing A: This case is similar to the (de) confinement properties of a planar 
AdS-Schwarzschild black hole. For any finite value of the dark radiation constant 
C > ("finite temperature"), eq. f JTTT) is satisfied, and the system is in the 
deconfined phase. For vanishing dark radiation constant C — 0, however, eq. flTUj) 
holds and the Wilson loop shows confining behaviour. 

3. The most interesting situation occurs for negative A: Depending on the value of 
the scale factor a (t) at any given time, either (I70|) or flTTj) can be satisfied, leading 
to transitions between confined and deconfined phases as the universe evolves. 
Physically, this transition is due to the competition between two effects: The 
screening effects of the thermal (dark radiation) energy always aims at driving 
the system into a deconfined phase, while the influence of the background cosmo- 
logical evolution on the static quark-antiquark potential can be either confining 
(for A < 0) or deconfining (for A > 0). Thus for negative boundary cosmological 
constant thermal screening and background-induced confinement can compete 
and result into (de) confinement transitions. 

For definiteness, let us consider the solution of (|T4|) for k = — 1 and A < 0, 




(70) 




(71) 




(72) 
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This solution thus describes an oscillating universe. Analysing (!70|) . ( ITT]) , we find 
two different regimes: Since the scale factor is bounded from above by a^max = 
1/y|A[, a large enough dark radiation constant C > R 2 /A always satisfies ( 17T|) . 
and hence quarks are always deconfined in this case. For C < R 2 /4:, however, the 
system oscillates between a deconfined phase at smaller scale factors ao(t), and 
a confined phase at larger values. For the marginal value C = R 2 /4:, the Wilson 
loop shows confining behaviour only at maximal extension of the universe. Our 
holographic setup thus describes a (supersymmetric) plasma with the qualitative 
properties observed in the evolution of our universe: Near the big bang singularity 
ao = the matter is in a deconfined state, and undergoes a confinement phase 
transition as the universe cools down. Figure [5] summarises the situation. 

These results, in particular in the latter case, need to be compared to the re- 
sults of [23], where it was argued that the Wilson loop is not a good measure for 
(de) confinement in AdS space. The arguments of [23J involve a conformal trans- 
formation between AdS space and half of the Einstein static universe (ESU), 
relating long distance behaviour in AdS space to the (universal) short distance 
behaviour in the ESU measured in turn by the Wilson loop. This argument fails 
in the backgrounds considered here since the conformal symmetry of Af = 4 SYM 
theory is broken by the gluon condensate, as well as by the conformal anomaly 
for A ^ 00 The results of sec. show that the Wilson loop is sensitive to both, 
and hence depends on the chosen conformal frame, measuring unambiguously 
the deconfinement properties of the chosen field theory state by coupling to the 
full energy-momentum tensor and to the gluon condensatejf) Similarly, we do 
not expect the arguments of [22] for a weakly coupled meron gas disordering the 
Wilson to simply carry over to strong coupling. 

The conformal anomaly is also responsible for the time-dependent nature of the 
holographic backgrounds presented here: The diffeomorphism relating these back- 
grounds to topological black holes [UJ [18] induces a conformal transformation on 
the boundary, which due to the conformal anomaly is not a quantum symmetry 
of the dual field theory. Thus, except for the case of vanishing boundary cosmo- 
logical constant, observables calculated in the two different conformal frames will 
generically be different, and therefore have to be calculated in the time-dependent 
background itself. We cannot resort to the equivalence with topological black 
holes to define a thermodynamic ensemble or to consider phase transitions, as 
e.g. studied for compactifications of CFT's on dS space-times |43j. Our setup 
hence is time- dependent and describes genuine non-equilibrium physics. We can 
however characterize the properties of the non-equilibrium state by calculating 
observables such as the Wilson loop via the machinery of gauge-gravity duality. 

For vanishing boundary cosmological constant the conformal anomaly vanishes, 
and the phase structure of the AdS-Schwarzschild black hole with flat (k = 0) 
horizon, which has 0(N 2 ) entropy density at finite temperature and 0(1) at 

7 This point was noted before in [42] . 

8 The gluon condensate is essential for confinement in the case A = 0, k = — 1, since the Wilson 
loop in the purely hyperbolic AdS-Schwarzschild black hole is screened [19] for all temperatures. 
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absolute zero [H] , coincides exactly with our results for the Wilson loop. For 
the hyperbolic (k = — 1) black hole, on the other hand, these two measures of 
confinement do not agree: From [T7J [18] it is clear that the black hole mass 
is given by // = C. The hyperbolic black hole has a nondegenerate horizon at 
r + = Lacis f° r A 4 — 0, and both its free energy and entropy are 0(N 2 ) at this 
point [IS]. Our Wilson loop on the other hand is confined at /i = C = 0. This is 
one of the rare examples where the density of states and the Wilson loop do not 
agree as measures of confinement, which can be traced back in this case to the 
effect of the gluon condensate which enters the Wilson loop via the string frame 
metric. 

As discussed in the introduction, different measures of confinement leading to 
different answers in curved backgrounds are not uncommon, and since the field 
theory is known exactly^ this case appears to be a good playground for a future 
investigation of the interplay of different measures of confinement We plan to 
come back to this point in a future work. 

Another interesting observation concerns the relation of the Wilson loop with the 
temperature (jlOl : Although this temperature has been derived in a adiabatic 
approximation, assuming the cosmological constant to be sufficiently small, the 
Wilson loop feels exactly this temperature, without any approximation. The 
reason is that the Wilson loop calculation in this section is done "locally in time" , 
i.e. by considering a string stretching into the fifth dimension at each fixed value 
in time, testing the presence of the horizon with (j40p . If the horizon is present 
the Wilson loop exhibits perimeter law, if not, the temperature (T4"U1) is zero or 
ill-defined, and the Wilson loop exhibits an area law. Thus, although ( 1401) can 
only be considered as an approximation, it exactly reproduces the Wilson loop 
behaviour. 



6 Summary and Discussions 

We have investigated properties of strongly coupled M = 4 U(N) supersymmetric 
Yang-Mills theory in the presence of a gluon condensate on cosmological space-times of 
Friedmann- Robert son- Walker type by applying methods of gauge/gravity duality. The 
dual gravity solutions are obtained from a Liu-Tseytlin Ansatz of IIB supergravity by 
solving the effective five-dimensional Einstein equations through a metric Ansatz first 
employed in the setup of braneworld cosmologies [T5] . By introducing a single inte- 
gration constant, the so-called "dark radiation constant" or "mirage energy density", 
Einstein's equations reduce to a single constraint equation. We solved this equation by 

9 The field theory is Af = 4 SYM theory with a gluon condensate [3]. This breaks conformal 
symmetry spontaneously, but not explicitly at the level of symmetry generators or Green functions, 
since the beta function still vanishes (see also |37j). 

10 This is particularly interesting in view of claims that the so-called 'precursor' states [3B], which 
create the 0(N 2 ) ground state entropy of the extremal hyperbolic black hole, are potentially relevant 
to the thermal screening of the quark- antiquark potential in the absence of the gluon condensate [19) . 
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Fig. 5: The (de) confinement transition as seen in the Wilson loop behaviour in an 
oscillating, open, universe (A < 0, k = — 1). During the periods depicted in blue, the 
Wilson loop shows confining behaviour, while during the red periods (close to the big 
bang singularity), the plasma is in a deconfined phase. The transition line, shown in 
black, is fixed by the value of the dark radiation constant C. 



introducing a boundary cosmo logical constant which a priori can be time- dependent. 
In this way the constraint is separated into a standard Friedmann equation and an 
equation for the sole undetermined function in the bulk metric. The resulting holo- 
graphic background is dual to M = 4 SYM theory in a FRW-type metric, and has a 
the gluon condensate. This bulk geometry determines the vacuum expectation value 
of the stress tensor. 

Holographically, the dark radiation term [201 El HH] induces a relativistic radiation 
component (of Stefan-Boltzmann form) in the stress-energy tensor. If the universe 
is static, this radiation component has the correct T 4 behaviour and N 2 scaling in 
the large N limit. If the universe is nonstatic, it is modified by a factor l/a 4 (t), as 
expected for a gas of relativistic particles. For nonvanishing boundary cosmological 
constant there is also the familiar conformal anomaly contribution proportional to the 
square of the boundary cosmological constant A 2 . 

Using the holographic stress-energy tensor, we have also clarified the holographic 
interpretation of a possible time-dependence of the cosmological constant A(i): The 
holographic stress-energy tensor is conserved if and only if the boundary cosmological 
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constant \(t) is time independent, A = 0. Requiring a conserved boundary stress- 
energy tensor is necessary to ensure that the dual field theory is boundary diffeomor- 
phism invariant. In general relativist's terms, the coupling between field theory and 
curved background geometry does not spoil the equivalence principle. Hence, with the 
standard AdS/CFT Dirichlet boundary conditions, only a time-independent boundary 
cosmological constant X(t) = A is holographically meaningful and consistent. 

The main result of this paper is the behaviour of the Wilson loop in the cosmo- 
logical background geometries, which we take as the measure of the (de) confinement 
properties of quark-antiquark pairs. We find an interesting interplay between dark 
radiation and cosmological constant: The dark radiation component drives the Wilson 
loop to deconfining behaviour, which may be understood as thermal screening of the 
quark-antiquark interaction. On the other hand, the boundary cosmological constant 
A can have deconfining or confining effect: Positive A (i.e. de Sitter like expanding 
cosmologies) drives the system into deconfinement, while negative A (i.e. anti-de Sitter 
like cosmologies) drives it towards confinement. For negative A there exist periodi- 
cally oscillating cosmologies in which the system undergoes periodic (de) confinement 
transitions: For dark radiation constants below a critical value set by the bulk AdS 
radius, the Wilson loop is deconfined when the universe is small (i.e. near the big 
bang singularity). After sufficient expansion it undergoes a transition to confining 
behaviour. This is in qualitative agreement with the expected behaviour in nature: 
Close to the big bang, i.e. at large temperatures, QCD matter should have been in 
a quark-gluon plasma state, and undergoes a confinement phase transition once the 
universe sufficiently expanded and cooled down. 

It should be stressed that in contrast to previous works [201 CEO HB], the standard 
holographic Dirichlet boundary conditions employed in this paper forbid bulk-boundary 
energy exchange. The dual field theory lives on a curved but fixed background geom- 
etry with no propagating graviton in the boundary theory. The Friedmann equation 
is obtained from the bulk Einstein equations. It is then natural that boundary diffeo- 
morphism invariance restricts the choice of bulk geometry. 

As a follow-up and in view of the interesting (de) confinement properties of the 
system described, it appears to be worthwhile to investigate the dynamics of fun- 
damental degrees of freedom. Such degrees of freedom can e.g. be introduced into 
these cosmological backgrounds by probe D7 branes [ID]. This would for instance al- 
low an investigation of the chiral symmetry breaking and its relation to Wilson loop 
(de) confinement, as well as possible chiral symmetry enhancement or suppression. 
Applying holographic renormalisation to these backgrounds will further clarify the re- 
lation between the free energy, entropy, and the properties of the Wilson loop, and thus 
shed light on the relation between these different measures of confinement. In a similar 
way the existence of a mass gap can be inferred from a fluctuation analysis around the 
bulk geometry. We will come back to these questions in a future work. 
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